The p-wave polaron 
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We consider the properties of a single impurity immersed in a Fermi sea close to an interspecies 
p-wave Feshbach resonance. We calculate its dispersion and spectral response to a radiofrequency 
pulse. In the presence of a magnetic field, dipolar interactions split the resonance and lead to 
the appearance of two novel features with respect to the s-wave case: a third polaron branch in 
the excitation spectrum, in addition to the usual attractive and repulsive ones; and an anisotropic 
dispersion of the impurity characterized by different effective masses perpendicular and parallel to 
the magnetic field. The anisotropy can be tuned as a function of the field strength and the two 
effective masses may have opposite signs, or become smaller than the bare mass. 



Understanding the physics of a single impurity in a 
degenerate ultracold gas has been essential in discover- 
ing the phase diagram of spin-imbalanced Fermi mixtures 
[1-3]. Close to a Feshbach resonance, the impurity state 
becomes truly many-body in character. A large effort has 
been devoted to its understanding, resulting in an impres- 
sive agreement between theory and experiment in recent 
years [4]. The emergent picture is that the quasiparti- 
cle formed by the impurity interacting with the back- 
ground gas can be fermionic (a "polaron" ) or bosonic (a 
"molecule" ) . 

So far, the polaron has not been studied close to higher 
partial wave resonances. Particularly interesting are p- 
wave resonances, as p-wave coupled supernuids are pre- 
dicted to display a richer phase diagram than their s-wave 
counterparts [5-8]. For example, the Bardeen-Cooper- 
Schrieffcr (BCS) and Bosc-Einstein condensation (BEC) 
regimes of these supernuids are separate phases [9] , and 
each of these can be either chiral or polar [7] . When con- 
fined to two dimensions, the superfiuid BCS phase is even 
topologically nontrivial [10]. Therefore it becomes cru- 
cial to understand the nature of the quasiparticles which 
could form such states. For 6 Li and 40 K, p-wave reso- 
nances between atoms in the same hyperfine state were 
found to have an extremely narrow magnetic width, of 
order < 1G [11, 12]. However, p-wave resonances in Li- 
K mixtures [13] proved to have larger magnetic widths, 
in the range 1 — 10G, indicating a stronger open-channel 
character of the p-wave interaction. Given the stabil- 
ity of magnetic fields in state of the art experiments 
(rj ±lmG), such resonances are now finally accessible 
for detailed study. 

Here we study the p-wave polaron problem by consid- 
ering a single impurity atom, labelled 4-, immersed in a 
spin-f Fermi sea, assuming that the two atomic species 
have the same mass and are strongly interacting in the 
p-wave channel, as may be achieved close to a p-wave 
Feshbach resonance. We will neglect any background s- 
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FIG. 1. (color online). Spectral function of the p-wave po- 
laron. The horizontal and vertical axes are — (k F v±i)~ 1 and 
ui/cf- The mi — ±1 and mi = resonances are located re- 
spectively at x = and x — S, with 8 =0(a), 0.5(b), 1(c), 2(d). 
Thick lines are the dressed molecules with mi — (solid) and 
m; = ±1 (dashed). Here — ko/kp = 10. 



wave interactions (see below). 

An interesting feature of p-wave Feshbach resonances 
is that they are generally split into a doublet according to 
the projection mi of the relative angular momentum onto 
the magnetic field axis. As first discussed in Ref. [12] , the 
mi = ±1 resonances are shifted towards higher energy 
from the mi = by the magnetic dipole-dipole interac- 
tion between the outer shell electrons in the presence of 
the external magnetic field. As we shall see, an important 
consequence is that when the energy splitting between 



the two molecular levels is larger than the width of the 
mi — molecule-hole continuum, a new polaron branch 
appears lying between the attractive and repulsive po- 
laron branches which had been previously observed close 
to s-wave Feshbach resonances [14, 15]. This is reflected 
in the spectral function of the impurity atom displayed in 
Fig. 1. Additionally, we find a regime close to the mi = 
resonance in which the effective mass of the attractive po- 
laron becomes negative (positive) along (perpendicular 
to) the magnetic field axis with the opposite behaviour 
near the mi — ±1 resonance. Clearly, these two new 
features could lead to many-body behaviour (e.g. collec- 
tive modes, possible superfluid states) which would be 
dramatically different from that of the s-wave system. 

The impurity may also bind to a particle from the 
Fermi sea to form a diatomic molecule, and as the in- 
teractions are increased this molecular state becomes en- 
ergetically favorable. We calculate the energies of the 
polaron and the molecule, and show that the position of 
the polaron-molecule transition shifts towards the BCS 
side for increasingly narrow resonances. 

To model the p-wave Feshbach resonance we use the 
following two-channel Hamiltonian [7, 16], and work in 
units where ft — 1: 



H = 
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For convenience we define bosonic operators b ll=x ^ y ^ z in 
terms of the closed channel I = 1 molecule operators 
&m,=o,±i as b x = (b + i+b-i)/V2, b y = -i(b + i-b-i)/V2, 
b z = b . We specialize to the case where the m; = ±1 
resonances are degenerate [17]. a a and a\ are the cre- 
ation and annihilation operators of fcrmions with mass 
m, and V is the system volume. The resonance splitting 
caused by dipolar anisotropy may be modelled by a posi- 



tive shift of the mi = ±1 molecule energies: e Xj . 
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Interactions couple the closed channel molecule to a 
pair of atoms in the open channel, and are described by 
the momentum-dependent coupling constant .g(|p|)p^- 
The coupling vanishes above a cutoff A of the order of 
1/R e , the inverse van der Waals length. In the low en- 
ergy limit, the physical results should not depend on the 
actual shape of g and we choose the cutoff function to be 
proportional to a step function, ,g(|p|) = .gO(A — p). 

In the absence of the Fermi sea, the bare propagator 
of the closed channel molecule can be read off from the 
Hamiltonian as 

£#,(p,«) = %77 — ^ rfffohO^,. (2) 

p ui — e M — p z /Am + lO p 

The renormalized molecular propagator D is dressed 
by the polarization bubble II as shown in Fig. 2a. In 




FIG. 2. (a) The renormalized propagator of molecules (thick 
wavy line). The thin wavy line is the bare molecule while the 
straight lines are fermions. (b) The diagram which leads to 
the T-matrix, Eq. (4). 



vacuum II" = II ^,, is diagonal and takes the form 



n°( P , 



dh 



q 9 (q) 



3 J (2ir) 3 u> — q 2 /m — p 2 /4m + iO 



(3) 



Thus the propagator in vacuum is also diagonal, i.e. 

DO I/ (p, w ) = W{[^(p^)]- 1 -n°(p, w )}. 

The elastic scattering between a spin-f and a spin- 4- 
atom in the ladder approximation is described by a T- 
matrix. The vacuum T-matrix in the present problem is 
given by (see Fig. 2b) 



v • 9 2 (k)kX, 

^ [D^(0, P/m)]- 1 - n°(0, fc 2 /m) ' 



E 



(4) 



Note that if the detuning e M is independent of ji the in- 
teraction is isotropic and the T-matrix is proportional to 
k • k'. However, the splitting of the p-wave resonance by 
the dipole-dipole interaction in the presence of a mag- 
netic field breaks rotational symmetry. In order to relate 
the parameters of the model to the physical obscrvables, 
we compare with the low energy expansion of the p-wave 
scattering amplitude /(k, k') = — (m/47r)T(k, k'): 
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(5) 



Here u M is the state dependent scattering volume and ko, 
with dimension of momentum, is the p-wave analogue 
of the effective range [18]. Evaluating the polarization 
bubble and inserting in Eq. (4) yields the relationships 
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with ci = m 2 g 2 A/6ir 2 . The vacuum molecule propagator 
is then 



dUp^) 
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—Vfi 1 + \k Q {muj - p 2 /4) - (p 2 /4 - muj - iO) 3 / 2 
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Note that the parameter kg is naturally large and nega- 
tive, of the order of the cut-off A [19, 20]. 

The p-wave resonances in the many-body system are 
characterized by two dimcnsionless parameters. Of these 



7 ~ (1 + C2)fci?/fco controls interactions at the scale of 
the Fermi momentum kp (i.e. many-body physics). Res- 
onances with 7 < 1 (7 » 1) are termed narrow (wide) 
and quantum fluctuations are suppressed for 7 <C 1 [7]. 
The second parameter c 2 controls interactions at the scale 
A (i.e. few-body physics) [19] and distinguishes strongly 
coupled (c 2 3> 1) from weakly coupled (c 2 <C 1) systems. 
Identical fermions interacting close to a strong Feshbach 
resonance may form trimer states [20]. However, using 
the method of Ref. [19], we have checked that trimcrs 
consisting of two identical fermions and an impurity of 
equal mass do not form in vacuum in the absence of res- 
onance splitting. 

We now turn to the question of the many-body state 
of a spin- J, impurity immersed in a spin-f Fermi sea at 
zero temperature. In the s-wave case both the attractive 
and repulsive polaron branches are accurately described 
by dressing the impurity with a single particle-hole ex- 
citation [2, 21]; we will adopt the same approximation 
for the p-wave case. Indeed, this approximation becomes 
exact far away from resonance, or in the limit of narrow 
resonances, 7 C 1 [7] . The propagator of the impurity 
with momentum p and energy to in the medium is 



Gi(p,u>) 



1 

u - p 2 /2m - S(p, oj) +i0' 



(8) 



with £ the self-energy. The energy of the polaron satisfies 
E = Rc[£(p,_E)]. The self energy in the single particle- 
hole approximation is given by the diagram in Fig. 3: 



E(p,w) 



d 3 c 



rn Fr (q)T(p,w;q,e q ) 



(9) 



where the Fermi function np^-(q) takes the value 1 if the 
state with momentum q is occupied, otherwise. The 
off-shell two-particle scattering T-matrix in the medium 






q\ (p q 



T(p, W ;q,e q )= 5 2 (|p-q|/2 // x . ; , , , ; 

x{[D cI (p + q,u; + e q )]" 1 -n(p + q,a; + c q )j \ (10) 

Here (p,w) [(q, e q )] are the f [l] atom momentum and 
energy entering the T-matrix, with e q = q 2 /2m. The 
polarization bubble in the medium is given by 



r d 6 ) 

n„( q ,c) = y — 



rf 3 k k^k v g 2 {k) [1 - n FT (k + q/2)] 



) 3 uj — q 2 /Am — k 2 /m + iO 

(11) 
Since |q| < A it is useful to write n = n° + (n - n°). 
Then, upon proper renormalization, the coupling in Eqs. 
(10) and (11) reduces to the bare value. 

The polarization bubble n is a tensor, as collisions with 
particles in the Fermi sea generally do not preserve the 
projection of the molecule's angular momentum. How- 
ever, D cl is diagonal, and the matrix inverse in Eq. (10) 
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p,u /' " p,w 

T„„(p,w;q,e q ) 



FIG. 3. The polaron self energy. 



may be performed as follows. We write the polarization 
bubble as 

n^(p + q,w + e q )=A^ + B(p + q) Al (p + q) y /|p + q| 2 . 

(12) 
In the matrix inverse of Eq. (10), sums containing only 
A or only B terms form geometric series. Resumming 
first the A terms, and successively adding the B terms, 
one finds 



^(p,o;;q,e q ) = 



t>. 



V. 



\p + <l\ 2 /B-V + 



, (13) 
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and 2? x = 



with T>± = "^2 
v^ (p+q)n(p-q)M 

The above formalism allows us to calculate the spectral 
function of the impurity atoms, defined as j4^(p,w) = 
— 2Im[(jj.(p, u>)}, which gives the spectral response of the 
impurity to a radiofrequency pulse of frequency u. In 
particular, the spectral function peaks at the energy of 
the quasiparticle states with a finite wavefunction over- 
lap with the bare impurity, i.e. the polarons. Our results 
for the spectral function close to resonance are shown in 
Fig. 1 for several values of the dimensionless resonance 
splitting S = 5 o12tt / (mg 2 k F ) . In the absence of reso- 
nance splitting, the picture is qualitatively the same as 
in the s-wave polaron case [2, 22, 23] with an attrac- 
tive (repulsive) polaron of energy lower (higher) than 
the bare impurity. However, an additional branch ap- 
pears for finite resonance splitting. For small resonance 
splitting, Sq < ej?, the intermediate branch disappears in 
the mi = molecule-hole continuum. For larger reso- 
nance splitting the intermediate branch shows up as an 
isolated spectral line which is a repulsive polaron close 
to the mi — resonance, approximately a free impurity 
between the resonances, and an attractive polaron close 
to the mi = ±1 resonance. 

In addition to energy, the polaron is characterized 
by its quasiparticle weight, the residue, given by Z = 
[1 — <9Rc{£(0, E)} /BE] evaluated at the quasiparti- 
cle energy. In Fig. 4 we display the energy and residue 
across resonance for various values of fco/fci?- We see that 
the absolute values of the polaron energies decrease for 
increasing \k n /kp\ as expected since particle-hole fluctua- 
tions become suppressed [7]. Indeed, in this limit the self 
energy (9), in the absence of resonance splitting (v^ = v), 
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FIG. 4. (color online). Energy (main) and residue (inset) 
of the two p-wave polarons for S = through the crossover 
for range parameters: — fco/fcF =30(solid), lO(dashed). The 
dotted lines are the perturbative result, Eq. (14), and the 
thin lines are the molecule energies as given by the Thouless 
pole. 
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FIG. 5. (color online). Inverse effective mass of p- wave po- 
larons moving perpendicular (solid) and parallel (dashed) to 
the magnetic field. Here —ko/kp =10, and the two resonances 
are split by S — 6, i.e., the mi — ±1 (mi = 0) resonance lies at 
x — (x = 6). The thin lines are the analytic weak-coupling 
results rn/m* L » w 1 -+- v± i \\k F /iv, valid in the regime where 



takes the form 



2(0, .E) 



1-Kvn 



fdq 



-v- 1 + ±mk (E + e q /2) + iQ 

(14) 
As mentioned above, the impurity may also bind a 
particle from the majority Fermi sea to form a bosonic 
quasiparticle — a molecule. In the "Cooper pair" ap- 
proximation (i.e., with no particle- hole pairs), the en- 
ergy E mo \ of the molecular state is given by the pole of 
the T-matrix at p = and uj = E mo \ + ep. The latter 
approximation is an upper bound to the real energy of 
the molecule, which becomes exact for \ko/kp\ ^$> 1. On 
the BEC side, E mo i tends to E^ — tp, with E\, the energy 
of the molecule in vacuum. The molecule energy is in- 
cluded in Figs. 1 and 4. The latter illustrates that as the 
resonance becomes more narrow (for increasing —ko/kp) 
the transition from a polaronic to a molecular ground 
state takes place further towards the BCS limit, as is 
the case close to a narrow s-wave Feshbach resonance 
[14, 24-26]. Note, however, that the limit of narrow p- 
wave resonances corresponds to small densities, while the 
opposite is true for s-wave resonances. 

The presence of a magnetic field which splits the res- 
onances also breaks isotropy since it differentiates be- 
tween mi levels. One consequence is the appearance of 
a strongly anisotropic dispersion of the impurity close to 
either of the resonances. At small momenta, the disper- 
sion may be written as E(p) — E(0)+p 2 /2mf,+p 2 L /2m* ± , 
where p» and p± are the projections of p along the mag- 
netic field, and on the plane perpendicular to it. We 
compute the effective mass tensor, showing how, on the 
BEC side, its magnitude is generally smaller than the 
bare mass of the particle. This is allowed since, in this 
region, the polaron is no longer the ground state. More- 
over, close to the mi = resonance m/mt becomes neg- 



ative while to/to^ is still positive, the opposite behavior 
occurring in the vicinity of the mi = ±1 resonance. This 
is illustrated in Fig. 5. 

We now make a few final remarks. In the above, any 
background s-wave scattering between the two atomic 
species has been ignored. This is justified provided the 
background scattering is small, in particular for a p-wave 
resonance close to a zero-crossing of an s-wave resonance. 
At vanishing scattering energy the s-wave channel dom- 
inates. However, in the present system the low energy 
scale is set by tp. Then the total cross section for scat- 
tering of two atoms in the s-wave channel is a s — Aira 2 
while the p-wave cross section is a p = 12ir\f p (kp)\ 2 [27]. 
If |wfc^,| _1 ^> \k /kp\ then p-wave scattering is dominant 
for \vk F \ > |fci?a| while in the opposite case we require 
\koa\ < 1. Both requirements are achievable provided fco 
is not too large. 

Early experimental studies of p-wave Feshbach 
molecules showed that these were short-lived, with life- 
times in the range 2 to 20 ms [11, 28-30]. This was ex- 
plained in Refs. [19, 20, 31] as due to relaxation processes 
as well as possible recombination into trimers. However, 
radiofrequency spectroscopy has been able to study in de- 
tail the complete spectral response of metastable many- 
body states of strongly-interacting fermions with an in- 
trinsic lifetime of only 30ms [14] including their polaronic 
and molecular branches. Furthermore, a recent proposal 
to modify the p-wave resonant interaction by coupling 
to a long-range excited state with the use of an optical 
Feshbach resonance predicts the possibility of suppress- 
ing three-body recombination [32], potentially allowing 
for longer lifetimes. 

In this Letter, we studied the properties of the p-wave 
polaron. These are observable by radiofrequency spec- 
troscopy. In particular, the anisotropic dispersion rela- 



tion may be studied using angle resolved photo emission 
spectroscopy. Our results are also relevant to studies of 
resonantly enhanced atom-dimer scattering near the ap- 
pearance of confinement induced p-wave trimers [33, 34] 
in polarized gases, in which case the dimer acts as the 
impurity investigated here. 
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